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Abstract 

We investigate the cosmological perturbation of two-scalar field model during the reheat- 
ing phase after inflation. Using the exact solution of the perturbation in long-wavelength 
limit, which is expressed in terms of the background quantities, we analyze the behavior of 
the metric perturbation. The oscillating inflaton field gives rise to the parametric resonance 
of the massless scalar field and this leads to the amplification of the iso-curvature mode of 
the metric perturbations. 
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I. INTRODUCTION 



Recently, a theory of reheating after inflation is developed and importance of the os- 
cillating scalar fields is recognized The non-linear interaction between the scalar 
fields amplifies the fluctuations of the scalar fields by the effect of parametric resonance and 
the energy of the oscillating field is transfered to the fluctuation of the massless field by 
catastrophic particle production. A large amount of studies about the reheating process 
has been done 0. From the viewpoint of the structure formation of our universe, one of 
the most important question is whether the growth of fluctuations of the matter field dur- 
ing reheating affects the large-scale inhomogeneities in the universe or not. To answer this 
question, we must study the evolution of the metric and the matter fluctuation in general 
relativistic treatment. 

Concerning the works on the theory of cosmological perturbation, several authors studied 
a single scalar field model in the context of the reheating scenario and investigated a role of 
the coherent oscillating field on the evolution of large-scale structure in the universe 0,^,^. 
Hamazaki and Kodama analyze the metric perturbation in the model with two-component 
fluid and discuss the effect of parametric resonance p. They evaluate the curvature per- 
turbation ( (Bardeen parameter) by replacing the scalar flelds with the perfect fluids. The 
effect of parametric resonance is all included in the energy transfer term which should be 
given by phenomenological description. They conclude that ( is well-conserved during re- 
heating and the parametric resonant decay does not affect the scenario of the structure 
formation. Although their conclusion seems to be valid in the old version of the reheating 
scenario, which is dominated by the Born decay process, it is not clear that their analysis is 
appropriate when the non-linear evolution of the scalar flelds may play an important role. 
We should keep in mind that the parametric resonance does not occur in the perfect fluid 
system. 

In this letter, we analyze the cosmological perturbation of the scalar fleld model to clarify 
whether the dynamics of coherently oscillating scalar fleld affects the metric perturbation 
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or not. The model we treat here contains two scalar fields 0, x with the potential V^(0, x) = 
m^0^/2 + g'^cjP'x^ /2 where is the infiaton and x is the massless boson field. 



II. BASIC EQUATIONS 

We consider a homogeneous, isotropic and fiat cosmological model with scalar type metric 
perturbations. In longitudinal (Newtonian) gauge, the metric is given by 

ds^ = -(1 + 2^)dt^ + a^{t){l - 2^)5ijdx'dx^ . (1) 

$, \1/ correspond to the gauge-invariant potential, which characterize the physical mode 
of the metric perturbations. The Bardeen parameter which represents the spatial curva- 
ture perturbation on the uniform energy density slice in super-horizon scale is written by 

ig,[n]0,[i3| 



C = $-^($ + i/-^«i'); H^-. (2) 
ri CL 

For the later analyses, we use the following dimensionless variables: 

T = mt, h = —, ip^ = J- , cp2 = J- , ^ = ^ • (3) 

m V 3 nipi V 3 nipi Att \ m / 

In terms of these variables, the background Einstein equations are given by 

h^ = ^? + v'i + 2U{^u^2y, U=^vl + ^Xvlvl (4) 
dU 

^': + Ship[ + — = (2 = 1,2), (5) 
where (') = d/dr. We use the following gauge-invariant variable to describe the perturbation: 

= Sy^i + (6) 

where 6ipi is perturbation of the scalar field ipi. The variable Qi was first introduced by 
Mukhanov |T^ and all the other gauge-invariant quantities can be expressed by these vari- 
ables. Using the fact $ = \E' in the present model, the Bardeen parameter defined by (^ is 
written by 



C = -72— 7^2 (^'l^^l + <^2Q2) 
V'l "T 7^2 



(7) 



The evolution equation for Qi is simply expressed by [U 

2 2 r 52f; « /.3 



A; 



ma 



6 / a" 



Q, = 0. 



(8) 



This is our basic equation. 



III. SOLUTION IN LONG WAVELENGTH LIMIT 

As we are interested in the formation of the large-scale structure in our universe, it 
is important to investigate the metric fluctuation outside the scale of the Hubble horizon. 
Fortunately, we can obtain the exact solution of eq.® in the long wavelength limit {k 0) 
and the result is expressed in terms of the background quantities. 

Let us suppose the solution of the background field is written as ipi{a, C), where a = log a 
and C is the parameter given by a suitable choice of the integration constants, a plays a role 
of the time parameter and we can trace the time evolution of the trajectory of the background 
solution in configuration space ip2) using this parameter. C distinguishes trajectories in 
this space and remains constant on a each trajectory. Using the time parameter a, the field 
equation (j^) is rewritten as 

+ (s + h-^—\— + h-^— = (9) 
da^ \ da ) da d^pi 

Differentiating with respect to a, C, we can show that the equation of motion for the 
tangent vectors di^i/da, dipi/dC becomes (||) in the limit k 0. Therefore these vectors 
are the independent solutions of in long wavelength limit [|19 . 



Since (H) are the second order coupled differential equations, we have four independent 
solutions. The remaining two solutions can be obtained as follows: define the component of 
a 2 X 2 matrix X as 



Each row vectors in X satisfy (|]) . Expressing the remaining solutions by the matrix Y and 
putting Y = X ■ P yields 

P" + (3/i + 2X-^X')P' = O. (11) 
We can obtain P' = X^^{X^^)'^ /a^ by using the fact that X' -X^^ is the symmetric matrix 



19| . Thus the general solution of (H) in long wavelength limit is given by 

dr 



2 



X 



-1\T 

where (^) represents transpose of the matrix and cf{j = 1,2) are arbitrary constants. 



(12) 



IV. BACKGROUND DYNAMICS IN REHEATING PHASE 

We can know the behavior of the long wavelength perturbation from (|T^) if we obtainthe 
background solution. As we pay our attention to the reheating phase of the inflaton domi- 
nated universe, we assume the condition ~ 1 and Ay^l ~ 1- The former is the condition 
for the oscillation of the scalar fields. The latter condition implies that the potential energy 
is dominated by the massive inflaton. 

We first solve the background equations in a naive treatment. Because Ayjg ~ 1? the 
evolution equations (^(^ can be solved separately and we get oscillatory behavior of the 
inflaton ipi ~ (^(r)cosr, where oc 1/r. Substituting this into (^ again, we find that the 
equation for the field ip2 becomes the same form as the Mathieu equation: 

(a^/V2)" + [A + 2q cos (2r)] (a^/V2) = ; a oc r2/^ (13) 

where the coefficients A and q are time dependent functions given hj A = 2q = A(^^(r)/2. 
The stability/instability chart of the Mathieu equation says that has unstable behavior 
atl/3<g<l, in the region of the first resonance band of the Mathieu function [|r^. 

To investigate the metric perturbation during the reheating, we must study the back- 
ground evolution more precisely. We shall analyze the background dynamics by using the 
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renormalization group (RG) method. RG method is a technique of asymptotic analysis which 
improves the resuh of the naive perturbation and it provides an unified approach to solve 



differential equations including singular perturbation method |T^ ||f 7|| . 

To apply the RG method, we identify the small expansion parameter e as the amplitude 
of the scalar fields (fi ~ 0{e) < 1. Since we want to take into account of the non-linear 
interaction, we assume A = eA ~ which implies A is not small parameter. Then the 
background quantities are expanded as follows: 

(pi = e(^i(i) + eVi(2) H , h = eh^i) + e^h(2) H • (14) 

The crucial treatment in the present model is to insert the mass term uj'^(p2 in the equation 
of (H)- We put uo"^ = 1 + ea and take the limit ea ^ — 1 after the end of calculation. This 
treatment is necessary to extract the effect of parametric resonance in field. 

In appendix, we apply the RG method to the background equations (§)(|^). The final 
form of the solution up to 0{e) becomes 

~ /icos 9 cos (r + 6*), (y92 — ^ sin 6 cos (r + ■?/'), (15) 

where we replace e/i(i) with h. The variables Q,9,ip and h are time dependent and deter- 



mined by so-called RG equations (|A^)- (|A7|) which are the result of the renormalization of 
the initial amplitude. 

We further impose the condition ^ V^i (© ^ 1)) which makes our analysis easier. We 
will restrict our analysis in the case of ^ 1 hereafter. ( |A6|) becomes ^' ~ and we can 



set 6^ = without loss of generality. Eqs. (|A5|) (|X7D become 

6' ~ ■^g(r) sin7 ■ 9, 7' ~ — 1 + g(r)(2 + C0S7), (16) 
h' = -lh\ (17) 

where we define the time dependent parameter 

g(r) = \xe. (18) 



Since ([T7|) yields h oc 1/r, g(r) is the decreasing function of time. Defining the new variables 
u, V hj u = Q cos (7/2), f = 6 sin (7/2), ( [T6| ) is rewritten as follows: 



fu'] 






2 


V 



J 



(19) 



l-q 
3g- 1 

The eigenvalue of the matrix in the right hand side of (^) is real for 1/3 < g < 1 and 
B can grow during this epoch. We see that the function q corresponds to the coefficient 
q in the Mathieu equation ([13|) . Recall that we have ifi ~ hcosr, (p2 — /i© cos (r + 7/2) 



from ([ISD , we recognize that the amplitude of (p2 has the growing behavior as the result 
of the parametric resonance. Soon after q becomes smaller than 1/3, O reaches a constant 
value. Fig.l shows the evolution of the background fields obtained by numerical integration 
of The field (the solid line) is amplified during the short period denoted by the 

horizontal arrow, which corresponds to the resonance epoch 1/3 < g < 1. We can evaluate 
how much is amplified through the resonance epoch. From (0), we have 



6 ~ 60 exp 



dr'q sin 7 



(20) 



Using this expression, we get the upper bound 



e(r^) 



< 



exp [C: 



VA); 



1 



1 



.|l--^^0.07. 



(21) 



where we used the relation q{Tf) = 1/3, qi^Ti) = 1. The above expression is valid if the 
condition 6 <^ 1 is satisfied throughout the resonance epoch. 



V. EVOLUTION OF THE BARDEEN PARAMETER 

Behavior of the background solution obtained by RG method is translated to the metric 
perturbation ( using the long wavelength solution (|T2|). The reduced system (|T6]) has the two 
integration constants: The initial values of 7, 0. During the resonance epoch, 7 approaches 
a constant value which does not depend on the initial value. We can regard Qq (the initial 



value of B) as the parameter C described in Sec.|lTli That is, 0o is the parameter which 



distinguishes the trajectories in the reduced system (6,7). Substituting ([T5| ) into (|T2D and 
evaluating the integral adiabatically, we find that the dominant contribution to the Bardeen 
parameter ( comes from the solution Q'^'^ = dipi/da and 

Qr = ^~-/ie,c sine COST, g-« = ^~/ie,c cose cos (r + |), (22) 

where Q^c = exp [/ dr'qjl ■ sin 7]. The other independent solutions are recognized as the 
decaying modes |T^ . 

The solution Q'f' is referred to as the adiabatic growing mode [jl8|. Substituting Q'f' 
into (1^), the Bardeen parameter remains constant in time. The constancy of C arises in the 



hydrodynamical perturbation without the entropy fiuctuation |jT2|. For the solution Q^*", 
(0) gives 

^00^ ,^ ^ ^ sin r COST — sin (r + ^) cos (r + ^) . . 

C-'' ^ hQ^c sin e cos e ■ ^ \ 2} (23) 

cos"^ e sm r + sm e sm [t + ^ j 

We see that has the periodically sharp peaks around the zero points of (f'l. This behavior 
also appears in the single field case |^. Except for the short interval of the peaks, 
traces the amplitude of Qf". Therefore deviates from zero when Q^c grows due to the 
parametric resonance. Note that vanishes when the background field ip2 is zero (e = 0). 
Because the initial amplitude of the curvature perturbation is vanishingly small, we can 
identify QY° is the contribution of the iso-curvature mode of the perturbation. 

We confirm the behavior of by solving numerically. In Fig. 2, we show the 

time evolution of the iso-curvature mode. In contrast to the constancy of the adiabatic 
growing mode (the dashed line), we observe that (the solid line) is amplified during the 
resonance epoch when the amplitude of the background field grows (see Fig.l). After the 
resonance epoch, approaches zero due to the Hubble damping which can be seen by the 
analytic result (P^D . Using (^T]) (^) , we can estimate the growth of C^" except for the sharp 
peak: 



< ^i^-exp(2c.v^). (24) 
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VI. SUMMARY AND DISCUSSIONS 



In this letter, we have investigated the cosmological perturbation of the model with 
two scalar fields in reheating phase after infiation. Using the background quantities, long 
wavelength solutions of the perturbation are obtained. Applying the RG method to the 
background dynamics in the coherent oscillating stage, we found that the massless field 
gets the effect of the parametric resonance. This leads to the amplification of the iso- 
curvature metric perturbations. The result implies that the Bardeen parameter deviates 
from the initial amplitude in reheating phase even though the initial amplitude determined 
in the infiationary stage is dominated by the adiabatic mode. Therefore the analysis using 
the scalar field model is essential to study the metric perturbations when the non-linear 
dynamics of the scalar fields plays an important role. 

The most important observed quantity is the power spectrum Vc^{k). Consider the evo- 
lution of the metric fiuctuation C^{k) whose initial amplitude is given by the quantum fiuc- 
tuation of the scalar fields during inflation. As long as the wavelength of the fluctuation 
is outside the Hubble horizon, time evolution of C,{k) is well-approximated by the solution 
for the long wavelength mode. The power spectrum for the present length-scale can be 
evaluated from the amplitude of C,{k) when the fluctuations re-enter the Hubble horizon. 
Soon after the metric fluctuation C^{k) suffers the parametric ampliflcation during reheating, 
it approaches to the initial amplitude due to the Hubble damping. The horizon re-entry 
time becomes earlier as the wavelength of the fluctuation becomes shorter. Therefore, we 
can expect that the effect of parametric resonance on the power spectrum at the horizon 
re-entry appears at the small scale (large k). We will give a detail discussion about evolution 



of the power spectrum in the forthcoming paper 19 
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APPENDIX: RG METHOD APPLYING TO THE BACKGROUND SYSTEM 



In this appendix, we apply the RG method to the background system (^(^. Following 
the treatment described in Sec.J^, we substitute ( |14|) to and expand in powers of e: 



0{e'): + ¥^,(1) = 0, (2 = 1,2) 



V'l{2) + Vl{2) = -3/l(i)V3i(i) - AV52(1)</31(1), 

V'2{2) + V^2(2) = -3/;.(l)V52(l) - 0-V'2(l) - Av3i(i)V22(l)- 



We first solve the above equations order by order. The solutions up to 0{e'^) become 



V9i = e 



+ e(r - To) ^-\hi)Ao + (2|5o|% + i?oX) 



and 



V?2 = e 



Bo + e(r - To) |-^/^(i)5o + (2|Aop5o + AIB*) + '-aB^ 

2^ a2u J3t 



+e'-Af,Boe''^ + c.c. + 0{e' 



where c.c. means complex conjugate and is given by 



^1(1) =2^|Ao|2 + |5oP. 



(Al) 



Aq, Bo are the amplitudes at the initial time t = tq. The above solutions are correct within 
e(r — To) ^ 1 because the secular terms which grow as (r — tq) makes the perturbative 
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expansion break down. To improve the perturbation series, we introduce an arbitrary pa- 
rameter fi. Splitting r — To as t — fi -\- fi — Tq and absorbing the terms containing /i — Tq into 
the renormahzation constant Za, Zb defined by Aq = y4(yu)Z^(ro, /i), Bq = B{^)Zb{tq, ii). 
This procedure is always possible for each order of expansion and we have the /i dependent 
solutions: 



V?i = e 



A + e{T- fi) \-h^i)A + (2\B\^A + B'^A*) 



+e^-B^Ae'^^ + c.c. + 0{e^) 
8 



= e 



B + e{T- fi) + {2\A\^B + A^B*) + '-aB 



+e'-A'Be'''^ + c.c. + C(e^ 



(A2) 



(A3) 



Since n does not appear in the original equations, the solution should not depend on /x. This 
requires {d(pi/dfi)r = for arbitrary r and we obtain the RG equations: 



dA 

dB_ 

dfi 



+ (2\B\^A + B^A' 



+ '-aB + I- (2\A\'B + A^B' 



2 " 2 2 

By setting /i = r in ( |A2D (|A3D , the secular terms are eliminated and we obtain the regular 
perturbation. Rewriting the variables A, B by 

A = cos 0e^^ B = sin Qe'^, 

and using ( |A1|) , we get the solution (|T5D and the above RG equations reduce to 

A 



O' = e — h,,) sin 7 sin 29, 
16 ' 

A 



7' = -1 + e^/iji) cos 26(2 + cos 7), 
^' = e^/i2 cos7sin2e, 7 = 2(^ - ^), 



(1) 



-^2^1)' 



(A4) 

(A5) 

(A6) 
(A7) 



where we set ea = —1. 
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FIGURE CAPTIONS 



Fig.l The amplitude of the background fields (pi and is shown as a function of time r by 
solving (^(^ numerically. The broken and the solid lines represent the infiaton(y9i) 
and the massless field(v92) respectively. During the period denoted by the horizontal 
arrow, the field ip2 is amplified by the effect of parametric resonance. We used the 
initial values which correspond to the expression (|T3p at r = with the parameters 
A = 4000, h = 0.035, 9 = 0.01, 7 = 0.4, 9 = 0.0. 

Fig. 2 The Bardeen parameter C^^" is plotted as a function of time r by solving 

numerically. We calculate with the same initial condition as described in Fig.l. 

As for the perturbations, we set the initial values corresponding to the expression 
( p2D at r = 0. The figure shows that C"*" (the solid line) is amplified during the 
resonance epoch denoted by the horizontal arrow and approaches zero due to the 
Hubble damping. We also plotted the Bardeen parameter for the adiabatic growing 
mode (the dashed line) obtained by the substitution of Q";'^ into (|^). 
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Fig.l 
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